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ABSTRACT 

We exploit the detection of three distinct stellar subpopulations in the red giant branch 
of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to 
resolve directly the evolution with radius of the dark matter mass profile. We find 
that a cored dark matter halo provides a perfect fit to the data, being consistent with 
all three stellar populations well within 1-sigma, and for the first time we are able to 
put constraints on the core size of such a halo. With respect to previous work, we do 
not strengthen the statistical exclusion of a dark matter cusp in Fornax, but we find 
that Navarro- Frenk- White haloes would be required to have unrealistically large scale 
radii in order to be compatible with the data, hence low values of the concentration 
parameter. We are then forced to conclude that the Fornax dwarf Spheroidal sits within 
a dark matter halo having a constant density core, with a core size of ro = llo^kpc. 

Key words: galaxies: kinematics and dynamics - Local Group - galaxies: individual; 
Fornax dSph 



1 INTRODUCTION 

Dynamical modelling of the dwarf Spheroidals (dSphs) 
has now been tackled with a number of different tech- 
niques. At the heart of these efforts lies the importance 
of constraining the properties of the central dark mat- 
ter density profile, so as to advance our understanding of 
galaxy formation towards the low- ma ss end of the mass- 
spectrum. J eans equation analyses fe .g.lWalker et al. 2009b; 
Lokas 2Q Q9|) , phase-space analyses ([Wilkinson et al.l 12002: 
WJ12007: A morisco & Evansll201lh as well Schwarzschild 
modelling (jjardel fe Gebhardtl 12012 : iBreddels etaD l2012h 
have all now been used in this field, with varied degrees 
of success. Nonetheless, it is fair to say that the main step 
forward has not really come through brute force in the mod- 
elling techniques, but rather from the realization that the 
presence of multiple stellar populations in dSphs could be 

turned into a powerful strength^ 

After the first work of iBattaglia et~ai1 (|2008h . the 
method has now been amplified to overcome some of 
the deficiencies of Jea ns analyses ([Evans et al ] 120091 : 
lAmorisco fc Evansll2012al ) and to give the required robust- 
ness to the identification and disentanglement of the co- 
existing stellar subpopulations (Walker & Penarrubia 2011, 
WP11 in the following). Very recently, lAgnello &; Evansl 
(|2012l ) (AE12) have brought the argument into sharp focus. 
They show that the additional evidence coming from co- 
existing subpopulations is extremely simple in nature, and 
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can be fully understood in terms of the independent require- 
ments to the global energetics of the system enforced by 
the different stellar subpopulations. By using the projected 
virial theorem (PVT) it is possible to identify the dark mat- 
ter profiles that are compatible with the photometry and 
kinematics of the subpopulations and exclude those in which 
the subpopulations would not all be in equilibrium at the 
same time. 

New rich and precise datasets will certainly allow and 
justify full use of more refined dynamical techniques in 
the near future, especially exploiting the spatially resolved 
prope rties of the line of sight (LOS) velocity distribu- 
tion ([Amorisco &; Evan s 2012B IJardel fc G ebhardt l20ll 
JG12 in the following). On the other hand, while the detailed 
orbital distribution of the stars is still a major uncertainty 
in the modelling, it is helpful to work within a framework 
that is entirely free from any anisotropy-driven degeneracy, 
such as the PVT. 

All the analyses mentioned so far point towards the 
presence of a constant density core at the center of both 
studied dSphs, Sculptor and Fornax. Together with pre- 
vious evi dence pertaining to low surface brightness galax- 
ies (e.g. Ide Blok et all 2001). these results have triggered 
a number of both theoretical and numerical studies. The 
aim is of course to understand the complex baryonic pro- 
cesses responsible for transforming the initial dark matter 
p ~ r -1 cusp predicted by Cold Dark Matter (CDM ) mod- 
els ([Dubinski fe Carlbergfll99ll : iNavarro et al.lll99/t ) 
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into a core. A number of suggestions have been made, 
including multiple epochs of m ass loss and gas re- 
accre tion (iRead Sz Gilmorel 120051 ). clumpy baryonic in- 
fall (|Cole et al.l 1201 ll ) and intermittent and concen- 
trate d star formation and su bsequent supernova feed- 
back (|Pontzen &; Governatollioi^ ) . Some or all of these pro- 
cesses may be able to deposit sufficient energy in the dark 
matter halo to erase the central cu sp. This picture has 
bee n tested with prom i sing r esults bv lTevssier et al. 
and lBrooks fe Zolotovl (2012), who show that a more realis- 
tic treatment of star formation and feedback can reduce the 
tension between CDM and the detection of cores in luminous 
dSphs. 

In this Letter, we revisit t he dark matter dens i ty pro - 
file of the Fornax dSph, after lAmorisco fc Evansl (|2012<J > 
(A 12) found that a three-population division better describe 
the data with respect to a two-population division (WP11). 
Both previous analyses of the Fornax dSph (WP11, JG12) 
were able to exclude the presence of an Navarro- Frenk- White 
(NFW) cusp with high statistical significance, although the 
size of the constant density core has remained poorly con- 
strained. WP11 obtain a measurement of the average loga- 
rithmic slope of the mass profile F (see eqn. fl3)) at a radius 
that is comparable with the half-light radii of their two sub- 
populations. This only allows them to put a lower limit to 
the core size. Similarly, JG12 find that a constant density 
core is required, but orbit-based modelling of a single stellar 
population would require a more extended set of kinematic 
tracers in order to put an upper limit to its scale length. 
By exploiting all the three different stellar populations in 
Fornax, we are able to find such a constraint. 



2 PROBING THE MASS PROFILE 

A 12 present evidence for the coexistence of three distinct 
stellar subpopulations in the red giant branch of the For- 
nax dSphs. These subpopulations have different metallici- 
ties, spatial distributions and kinematics. We refer to their 
Table 1 for details, and record here for convenience the rel- 
evant half-light radii: 



->MP 

^ M = 610 ± 25pc, 

->MR 



935 ± 65pc, 
610±25pc, 
437 ± 55pc, 



(2) 



where MP, IM and MR stand, respectively, for metal-poor, 
intermediate-metallicity and metal-rich. Fig. 1 displays the 
LOS velocity dispersion profiles of the three stellar sub- 
populations. These have been obtained by using the max- 
i mum l ikelihood technique described in lAmorisco &; Evans] 
(|2012bh . which takes into account the probabilities of mem- 
bership of each star to the identified stellar subpopulations. 

With this information, we produce a total mass estimate 
for each subpopulation . We us e the mass estimator proposed 
bv lAmorisco fe Evansl (|2012aT ) 



M [(1.67 ± 0M)R h ] = (5.85 ± 0.2) 



G 



(3) 



which has been proved to be accurate on the wide family 
of models with Michie-King phace space distribution. We 
also notice that use of alternative mass estimators with the 
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Figure 1. Metallicity distribution (in the Magnesium index - in 
A) and line of sight velocity dispersion profiles in circular annuli 
for the three stellar subpopulations identified in Fornax. Color- 
coding corresponds to increasing metallicity. The profiles use the 
same number of tracers per bin, although this is different for each 
population (34 for MR, 283 for IM and 123 for MP). 

sam e functional structur e - as for example the one proposed 
bv IWalker et al.1 (|2009bh and used in WP11 - would give 
identical results for the logarithmic slope of the mass profile 
= log Mi — log Mj 

logfl h>i -logfl h>j ' 1 J 

where i and j identify any two stellar subpopulation. The 
left panel of Fig. 2 illustrates the one and two sigma contours 
for each of the total mass estimates, with the same color- 
coding as in Fig. 1; the right panel of the same Figure shows 
the probability distributions of the three logarithmic slopes 
obtained by combining in pairs the distinct subpopulations. 

For any realistic dark matter density profile, the log- 
arithmic slope r(r) is a monotonic function of radius and 
decreases to zero at large radii. For a cored halo, T < 3 
within the scale radius of the core. For a cusped dark mat- 
ter profile p ^ r -7 , T < 3 — 7 near the center. Therefore, for 
an NFW halo, Y is never higher than 2. WP11 have mea- 
sured a value of Y that is systematically higher than 2 near 
the center of Fornax, and were able to exclude the presence 
of a cusp (7 ^ 1) with high significance (s ^ 96%). Since 
we have three different mass estimates, we can measure the 
logaritmic slope of the mass profile in multiple locations, 
hence assess directly its evolution with radius. In the ra- 
dial interval defined by the IM and MP stellar populations, 
which is the farthest from the center, the slope is centered 
on T « 2 (blue-purple curve in the right panel of Fig. 2). 
However, as soon as the radius is reduced by considering the 
MR stellar population, the mean value of Y sistematically 
increases. By considering the MR-MP pair (red-blue curve) 
we obtain (r) ~ 2.4, while, near the center, the MR-IM pair 
(red-purple curve) suggests (r) « 2.65. This is probing di- 
rectly the radial variation of the mass profile of the Fornax 
dwarf. 
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Figure 2. Left panel: estimates of the total enclosed mass M(r) as obtained from each of the three stellar subpopulations; color-coding 
is the same as in Fig. 1. Right panel: probability distributions of the logarithmic slope T obtained by pairing the mass estimates in the 
left panel (colors are paired accordingly). 



In a purely statistical sense, the exclusion of an NFW 
cusp we obtain from this analysis is not as strong as the 
one obtained in WP11: the tails of the probability distribu- 
tions of both Tmr,mp and Tmr,im suggest that the proba- 
bility of a density profile as steep or steeper than an NFW 
is p < 22%. Nonetheless, it is worth noticing that such a 
cusped density profile would be unrealistic for different rea- 
sons. For instance, the dashed curve in the first panel of 
Fig. 2 represents the best fitting NFW halo with scale ra- 
dius ro = 3 kpc. The agreement with the mass estimates 
is marginal and it is in fact impossible to obtain a better 
agreement for any NFW density profile with a smaller scale 
radius (hence even smaller (r) at the same radi i ) . Different 
studi es of Fornax (jlrwin &; Hatzidimitrioul 1 19951 ; I Cole et all 
20li WP11) suggest that its tidal radius is likely to be 
smaller than 3 kpc, hence presenting an evident incompati- 
bility with the required scale radius. Furthermore, even an 
NFW halo with a scale radius of 'only' ro > 2 kpc would al- 
ready have a concentr ation c < 9, well be low the expectation 
of CDM models (e.g. iMaccio et al1l2007l ). As a comparison, 
the full line in the lef t panel of Fig. 1 represents a Burkert 
profile (|Burkertlll995h 



/>Bur(r) - 
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with a core-size of ro ~ 1.4kpc, which provides excellent 
agreement with the mass estimates, being consistent with 
all three stellar populations well within 1-sigma. 



3 PROJECTED VIRIAL THEOREM 

We can gain a more systematic insight by making use of the 
PVT 



2K los + Wic 







(6) 



where, in the usual notation, K\ os and W\ os are the projected 
components of the pressure and potential energy tensors - 
explicit for mule are recorded in AE12. The PVT provides a 
fundamental means of exploring the energetics of a system 
with multiple stellar subpopulations. As in AE12, we assume 
that the photometry of the different subpopulations is well 



represented by a Plummer profile; we have performed the 
same analysis using exponential profiles and find consistent 
results. 

In the ideal case of no observational uncertainties, use 
of the PVT proceeds as follows. For an assigned dark matter 
density profile, once the surface brightness of a stellar sub- 
population p(r) and its velocity dispersion profile a\ os (r) 
are given, the PVT yields the characteristic density po that 
needs to be coupled with the scale radius ro in order to have 
an equilibrium configuration. In our application, we use a 
Monte Carlo procedure to propagate any observational un- 
certainties in the half-light radius of the photometry from 
eqn. ([2]) as well as in the kinematics (see Fig. 1), and con- 
struct a complete probability distribution for the character- 
istic density po at any fixed characteristic radius ro- 

The left panel of Fig. 3 shows the one-sigma virial 
stripes that we obtain when each subpopulation in embed- 
ded in an NFW dark matter profile (eqn. (pQ)); the right panel 
illustrates the case of a Burkert halo (eqn (0). In the NFW 
case, the one-sigma regions related to the IM and MP sub- 
populations overlap before ro = 1 kpc. This is in agreement 
with the previous Section, since we found (Tmrjm) = 2.0 at 
large radii. However, at smaller radii the MR subpopulation 
does not fit into this picture. Any NFW halo that is compat- 
ible at the one-sigma level with all three subpopulations has 
ro > 2kpc, with the difficulties mentioned earlier. On the 
other hand, for a constant density core, the three one-sigma 
stripes show a consistent overlap region. 

Use of a single subpopulation does not allow any infer- 
ence on the scale radius of the core if taken by itself. It is 
by multiplying the probability distributions defined by the 
three virial stripes that we can constrain the core size. Fig. 4 
displays the one- and two-sigma confidence regions associ- 
ated with the joint likelihood, together with the marginal- 
ized probability distributions for the characteristic density 
po and core size ro . All full lines are associated with a Burk- 
ert dark matter profile (eqn. J5J)), but we find that the data 
do not show any significant preference for other functional 
prescriptions of the core. For example, we have explored the 
density profile used by AE12 

Po 
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Figure 4. Joint likelihood analisys of the virial stripes: one and 
two sigma contours in the halo plain (ro; po) and marginalized 
probability distributions for a Burkert density profile (full lines). 
The short-dashed probability distribution is associated with the 
cored profile of eqn. (7), while the long-dashed distribution is 
indicative of an NFW halo. 



and we find no difference in either the virial stripes or in the 
maximum value of the joint likelihood. In fact, the dotted 
line in the lower panel of Fig. 4 has been obtained by using 
the density profile J7} , but this is virtually indistinguishable 
from the one obtained for a Burkert profile. As a comparison, 
the dashed line in the same panel is associated with the joint 
likelihood of an NFW halo. Increasing the scale radius ro is 
the only way to increase the average value of T in the radial 
region sampled by the three stellar subpopulations, therefore 
the marginalized probability distribution is a monotonically 
increasing (respectively, decreasing) function of ro (of c) . We 
then conclude that the 68% confidence interval for the core 
size is ro = lto ^kpc. 



4 CONCLUSIONS 

The Fornax dwarf is the brightest, and probably the most 
massive of all the dwarf spheroidal (dSph) companions sur- 
rounding the Milky Way Galaxy. It has a complex evolution- 
ary pat hway, as is already evident from its star form ation 
history (jColeman fe de JondbdoSlde BoereT al. 2012). The 
detec tion of localized overdensities ([Coleman et al.l |2004 
2005) as well as a complex dependence on metallicity of 
its rotational properties (|Amorisco &; Evansll2012cr ) support 
the idea of a relatively recent interaction with a smaller 
dwarf. Large spectroscopic samples o f thousands of red g i- 
ant stars in Fornax are now available ( Walker et al . 2009a), 
and have recently been separated into three populations 
- metal-rich, intermediate- metallicity and metal-poor - by 
lAmorisco fe Evansl (l2012cT ) using chemo-dynamical meth- 
ods. 

In stellar dynamics, multiple populations provide a 
powerful method to constrain the common gravitational po- 
tential in which the stars move. Here, we have exploited the 



division of Fornax into three populations to show that the 
stars reside in a dark matter halo having a constant density 
core with size of ro = l^Q^kpc. This is in contrast to the 
case when just the single, combined stellar population of For- 
nax is studied. Then, the k inematic data are equally consis- 
tent with cusps and cores (Evan s et al ] 120091 : IStrigari et al. 
l2010h . Although often cited as an issue, the ellipticity of the 
luminous dSphs is unable to alter the co nclusions on their 
dark m atter density profiles, as shown bv lAgnello &; Evansl 
(|2012h . Even an oblate NFW halo would still need to be so 
large and diffuse to be substantially different from expecta- 
tions. 

All this adds weight to the evidence that the dSph 
galaxies of the Milky Way possess dark matter cores. This 
evidence includes the su rvival of kinemati cally cold sub- 
structure in Ursa Minor (|Klevna et al.ll2003h , the longevity 
of the globular cluster population in Fornax (iGoerdt et al.l 
2006; Sanchez -Salcedo et al. 2006: ICole et al. ]|2012h, as well 
as the kinematics of stars i n Sculptor feattaglia et aT1 l2008 ; 
lAmoris~ fe Evansl l2012al: lAgnello fe Evansl l2012h and in 
Fornax (Wal ker fe Penarrubiall201lh . There is no direct ev- 
idence in favour of a dark matter cusp in any dSph. 

The nature of the central parts of dark haloes of dSphs 
is important for two reasons. First, indirect detection ex- 
periments are monitoring the sky for signals of the ex- 
pected 7 ray emission produced by pair annihilation of 
weakly interacting dark matter particles. For example, the 
Large Area Telescope (LAT) on the Fermi satellite has con- 
strained the integral 7 ray flux from 1 4 dSphs to place 
limits on the annihilation cross-section (jAbdo et all [2010; 
lAckermann et aL I l201lh . The beam -size of the LAT is com- 
parable to the angular half-light radius of most dSphs (~ 
0.5°). In turn, this is comparable to the angular scale where 
uncertainties f rom the dark matter profile on the J facto r 
are minimized (| Walker et al . 2011; C harbonnier et aDl201lh . 
Therefore, whether cored or cusped models are used makes 
only modest differences to the interpretation of the Fermi 
LAT data. However, arrays of atmosphe ric Cerenkov tele- 
scopes like MAGIC and HESS (e.g., lAleksic et aTll201ll : 
lAbramowski et all 1201 lh provide superior angular resolu- 
tion, and now the existence of a dark matter core rather than 
a cusp substantia lly reduces the expected 7 ray flux (e.g., 
lEvans et alll2004h . 

Secondly, it is a unique prediction of cold dark mat- 
ter theories that the centres of dark haloes are cusped. The 
dSphs of the Milky Way are overwhelmingly dark matter 
dominated, and so provide natural testing grounds for as- 
sessing this hypothesis. The finding that the haloes of at 
least some of the dSphs are cored may be interpreted in one 
of two ways. Either it provides evidence for warm dark mat- 
ter particles which - judging from the ph ase space density in 
the F ornax core - have mass ~ 0.5 keV ()Dalcanton &; Hoganl 
l200lh . Or, it provides evidence that baryonic processes 
such as supernovae feedback or episodic mass loss or in- 
fall have modified the pristine dark matter cusp. There is 
no shortage of sug gestions of processes th a t could eradicate 
the cusps (see e.g . Read fc Gilmorel 120051 ; ICole et aT]|201ll ; 
IPontzen fe Govern ato 2012). The increasing size and rich- 
ness of datasets on the dSphs, combined with the increasing 
sophistication of analysis tools, gives hope to the suggestion 
that observational data may soon allow us to discriminate 
between some of these possibilities. 
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Figure 3. Left panel: virial stripes for the three subpopulations embedded in an NFW gravitational potential. Right panel: virial stripes 
for the three subpopulations embedded in a Burkert gravitational potential. In both panels stripes are associated with 1-cr regions and 
color coding is the same as in Fig. 1 and 2. 
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